Abstract 



The hydrodynamic lift force that polymers experience near boundaries is known 
to be a crucial clement when considering rheological flows of dilute polymer solutions. 
Here we develop theory to describe the hydrodynamic lift force on extended polymers 
flowing near flat surfaces. The lift force is shown to display a non-monotonic character 
increasing linearly with the distance to the wall Z in the near-surface regime defined 
as Z < L, with L being the contour length of the polymer. At heights Z ^ L the lift 
force displays a maximum, and for Z > L we recover the well known far-field result in 
which the force decays as Z^^. Our analytical theory has important implications in 
understanding adsorption, desorption, and depletion layers of highly extended objects 
in flow. 
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There is a " depletion zone" in flowing solutions of macromolecules which is characterized 
by a drastic reduction in the concentration of the polymer near a surface. [5l [151 [121 [2] Over 
the past decade it has become clear, after early work by Jendrejack, et al., that this 
behavior is due to a hydrodynamic lift force, which is a direct consequence of polymer 
stretching under the influence of an external fluid flow. [151 [HI [T2l \TU[ [9l [6l [11] The origin 
of this lift-force arises from a no-slip boundary condition at a solid surface that breaks the 
fluid flow symmetry around a chain under tension, and the chain feels a resulting upwards 
force. This force has been widely investigated and applied to a variety of contexts, in 
particular the description of the aforementioned depletion zone and the case of polymer 
adsorption. [ISl[Hl[Il[ini[ni[IZl[Zl[Hl[Iil[2ni[I21[2I]Maand Graham were the first to consider 
this behavior in the context of adsorption, and showed how a dumbell polymer model could 
be developed that displayed the qualitative characteristics seen in experiment. |15] This 
model was further developed by Hoda and Kumar who extended this model to incorporate 
the effects of electrostatic interactions. |9] Theoretical and simulation work by Sendner and 
Netz has further investigated this phenomenon in the context of more refined bead-spring 
models that more explicitly incorporate the correct hydrodynamics. [20 l 118 1 fl ^ I19| Despite 
these advances in the theory, comparison of many of these results with experiment reveal 
qualitative, but not quantitative agreement. Larson, et al. has indicated that the depletion 
zone near a surface is considerably smaller than predicted by theory, by as much as an order 
of magnitude. [5] 

In this Letter we provide a general theory that describes the hydrodynamic lift behavior 
of a polymer at an arbitrary distance from the surface. Our theory provides an explanation 
of the discrepancies seen between experiment, simulation, and theory for surface depletion, 
and also an understanding of the lift forces that are relevant for polymer adsorption and 
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desorption. In particular we note that previous models addressing this force, such as 
the dumbbell model, uses an explicit far-field assumption, the validity of which has been 
questioned by Hernandez-Ortiz, et al. and Hoda, et al.. However, a complete quantitative 
description of this regime has not yet been developed. [6l [9] Here, we provide analytical 
results that capture the fundamental physical picture in both the far-field and near-surface 
regimes. Our analytical theory supposes no underlying chain model, and relies only on 
geometric parameters and flow conditions. It shows for the first time that the lift force 
on an extended object, like an extended polymer, can become non-monotonic, increasing 
linearly with the distance from the surface. In the far-field regime we recover the well 
known results that the force decreases quadratically with the distance. Of further interest 
is the fact that the crossover from the near-surface to the far-field regime where the lift 
force is maximal occurs at a height ~ L, where L corresponds to the polymer extension 
length. Since L can be extremely large, particularly in strong flows, the near field regime 
may extend far away from the surface. 

To develop this model, we only need to consider the geometry of a polymer chain near 
a surface. Here we provide a description similar to the blob model of a polymer chain 
under tension by representing the chain as a series of hydrodynamic beads whose length 
is equal to the ratio of the average length to the average width of the elongated polymer 
chain 2N = (AX) /(AZ) (see Figure [T^ for a schematic) . [Ij This description of the polymer 
chain is amenable to the incorporation of molecular theory; for example a straightforward 
calculation of this geometry in shear flow based on the work of DeGennes is performed in 
the Supplemental Material (SM), and yields the relationship 2N = ^\ + {'^tzY / {'2E'^), 
where ^tz = Wi is the shear rate rendered dimensionless through the use of the chain 
relaxation time tz and E is Peterlin's chain stretching parameter that is 1 at low extension 
and diverges at high extensions. [3j We introduce a single length scale 2a = (AZ) to describe 
the chain dimension normal to the wall that effectively corresponds to the bead radius. This 
consequently allows us to represent our system as a bead-spring chain that lies parallel to 
the surface, in the same spirit as Sendner and Netz (see Figure [Tja for a schematic). |2Ul 
[T5] We can then determine an analytical form for the hydrodynamic lift force of this 
fully elongated polymer using the Blake-Oseen tensor, which has been shown to provide a 
good description of hydrodynamics near surfaces in the regime of interest {Z > a) . [1] IT3] 
This geometry conveniently enables us to use a simplified version of the hydrodynamic 
interaction tensor, and captures the essential physics of the lift force in the limit of high 
stretching. This approach is also the relevant case for a desorbing or adsorbing polymer, 
which exhibits extended conformations as it interacts with both the external fluid flow 
and the surface. [201 E] Understanding a polymer in such a limit is also important when 
considering the behavior of a polymer in channels that are characterized by small length 
scales. 

The lift force on a polymer, or any other deformable component dispersed in a liquid 
medium near a surface, is due to the hydrodynamic interactions between individual com- 
ponents and the no-slip boundary condition at the surface. The simplest manifestation of 
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Figure 1: (a) Simulation data showing the accumulation of bead positions over 4800 r 
at 7r = 2.1 for a 0-chain of 100 beads. We approximate this geometry using a bead- 
rod model, which groups portion of the chain into beads (orange circles). This geometry 
is shown in (b), and has 2N = {AX)/{AZ) beads of radius a = {AZ)/2 parallel to the 
surface at height Z. Individual dumbbell-pairs feel a net force Ff which is equivalent to the 
gradient of the tensile force along the chain in the continuum limit. All pairs of dumbbells 
contribute to an overall lift force Fl- 

the hydrodynamic interaction between an entity and its surroundings is given by the Oseen 
mobility tensor ^jj^o(i"i) rj), which describes the entity as a single point force. [UH] This 
mobility tensor describes the effect of a point force Fj (r j ) on the surrounding velocity field 
Vj(rj) through the equation Vj(rj) = ^ij^oi^i^'^'^j) ' Fj(rj). To account for the effect of a 
nearby surface, where a no-slip boundary condition must be maintained, Blake introduced 
an image system that accounted for the aforementioned boundary condition to produce a 
new mobility tensor See the SM for more discussion of the mobility tensors. 

In our model we consider a finite chain of 2N hydrodynamic beads at a distance Z 
from the surface, each of radius a. Using these parameters, we can render all variables 
dimensionless (rescaled variables denoted by tildes) in terms of distances a, energies kT, 
and times tq = Qirrjo? / {kT) , where t] is the solvent viscosity. This geometry is shown 
schematically in Figure [T|3. The tension force Ft along the chain is considered to only be 
in the x-direction, such that the chain is in equilibrium in this configuration. Since we 
are primarily concerned with the lift forces in the direction normal to the surface, we only 
need to consider the resulting velocity field in the z-direction. The relevant component of 
the mobility tensor is then just fJ-ij ^iT^x) = (-3ra;Z^)/(27r?7(r^. + 4Z^)^/^). To evaluate the 
overall force, one can just define the effect of a given dumbell centered at the mid point 
of the chain that is separated by a distance 4na on the velocity field at a dumbell that is 
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separated by a distance 4ma, where n and m are just indexing parameters (see Fig. 1). 
The resulting hft force from the n dumbbell on the m dumbbell can be writen as: 



Fdin,m,Ff, Z) 
[in 
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where i^j is the magnitude of the force on the nth bead and the rescaled height Z = Z/a. 
To determine the overall lift force due to a given dumbell on the entire chain -FL,d(n, Ff, Z 
(including the dumbell itself), we sum over all dumbells m = 1 to A^: 
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where we have used the approximation that N is large enough that we can replace the 
sum by a continuous integral from to N . To find the total lift force, we perform the 
summation over all dipoles n and again replace it with an integration: 
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Since Ff, the overall force on the dumbbell, is the difference of the tension on either side of 
the bead, we make the replacement of this force with its continuum analogue, [SFt / dn). 
This represents a "kernel" in which the form of the load profile is input into the theory. 
In this paper we consider a profile that corresponds to shear flows with the extended 
polymer at a small angle 9 (another example is given in the SM). We incorporate the 
result {OFx/dn) ~ 6-7rr/a^7sin (20)?i. Performing the integration, we get the result: 
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For the far-field, we retain the proper Z ^ scaling: -Fl, shear ~ §7''" sin (26*)^, however 
below a crossover height Z* N there is a near-surface regime where the lift force becomes 
non-monotonic and the full equation [4] must be used. This theory can now be related to 
any microscopic theory of choice, with the geometrical parameters (A^, a, r) corresponding 
to chain parameters. We show a straightforward example of this in the SM using Gaussian 
dumbbells to obtain the well-known far-field scaling relationship for the depletion length, 
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Figure 2: (a) Load profile that corresponds to a chain in shear flow, which is placed into 



X]?^ ^j,x/ (2A^)). (b)Graph of Fl versus Z for a shear load profile. Solid lines show theoret- 
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ical results, dashed lines show the effect of bead-discretization at low values of Z (see SM), 
and dotted lines demonstrate the far-field dumbbell results. Simulation data is also shown, 
with filled symbols representing data without fluctuations and open symbols representing 
data including fluctuations. These results collapse onto a single curve (inset) using the 
scaling Z ^ Z/iV and F ^ F/iV^^ 
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Or results are confirmed by Brownian Dynamics simulations (details in the SM) on 
chains that are extended parallel to a surface. These bead-spring models represent the 
effective blobs that make up an extended polymer chain. For these simulations, we in- 
corporate the appropriate hydrodynamics by using the Rotne-Prager-Blake tensor and 
consider the chain with both the absence and presence of thermal fluctuations. |16| |13j For 
these simulations we use the force loading proflle that corresponds to shear flow, which is 
given by fl{dFf{ri)/dn) = 7rsin26'(fi,x — X^j^ fj,x/(2A^)) (shown in Figureji^). 

In Figure we plot the result of equation |4] for a shear flow for chains of 2N = 20, 
50, and 80 with an effective shear rate of 7Tsin20 = 7Te// = 0.1. We plot two types of 
simulations; open symbols indicate simulations where fluctuations are turned off, and fllled 
symbols indicate simulations with thermal fluctuations included. It is clear that, especially 
in the limit of large N, the theory and simulation match well. This plot demonstrates a 
non-monotonous lift force with two characteristic regimes for this loading proflle; the near- 
surface regime demonstrates an increasing Fl with Z, and the far-fleld regime demonstrates 
a ~ Z~'^ decay that is expected from the far-field dumbbell model (plotted as dotted 
lines in Figure I2r,) . We can rescale this graph for all values of N by rescaling Z ^ Z /N 
and Fl — t- Fl/N^ to collapse all of the curves for a given value of 7Te// onto a single curve. 
This is shown in the inset of Figure [2] There are deviations between theory and simulation 
at small values of Z due to the discretized representation of the chain in the simulations as 
opposed to continuous representation of the chain in the theory. A discretized version of 
our theory can be developed that only requires a few terms at low Z (see SM) which more 
accurately reflects the low Z behavior in our simulations, and is shown as a dashed line in 
Figure ^p. We expect that a real chain would more closely resemble the continuous case, 
as the discretization of the chain is arbitrarily done to develop a convenient simulation 
model. 

The size and shape of the depletion layer near a surface can now be calculated using 
equation This is done by describing a potential of mean force Ul = -J^ FlcLZ. If the 
only force acting on the polymer in the Z-direction is the lift force, we can write the concen- 
tration c{Z) of a dilute polymer solution at height Z from the surface as c{Z) = c{oo)e~^^ . 
We plot sample proflles for our model, for the dumbbell model (where we use both the 
far-field result and the dumbbell result given by Hoda and Kumar) for shear fiow in Fig- 
ure [D] There is a striking difference between these two profiles, with drastic differences 
appearing due to the significantly lower lift forces close to the surface using our elongated 
chain model rather than the dumbbell model. We also predict a finite concentration at the 
wall, at least until steric or other short-range wall forces become significant, since the effec- 
tive lift force potential does not diverge at Z = 0. These characteristics have been noticed 
before in experiments, with depletion layer widths that are smaller than the predicted ones 
by as much as an order of magnitude, and which are often accompanied by finite surface 
concentrations of polymer. [5j Our model clearly provides an explanation of the origin of 
the difference between experiment and previous theories, which do not adequately account 
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Figure 3: (a) Concentration profiles highlighting the difference between the elongated chain 
model (solid lines), the dumbbell model (dashed lines), and the far-field dumbbell model 
(dotted lines) for a number of different shear rates near a surface with = 10 (given in 
terms of geometric parameters). The depletion region is often considerably smaller when 
near-surface hydrodynamics and elongated chain geometries are considered. While direct 
comparison of these results to Weissenberg numbers is dependent on the specific chain 
model, a simple scaling model (see SM) using these values suggests they correspond to 
Wi ~ 1 — 10. All three cases correspond to the same far-field behavior, (b) Probability 
distribution function for chains in a sheared slit flow. The near-surface hydrodynamics, 
which are considered in the elongated chain model, have the effect of greatly decreasing 
surface depletion in channel widths on the order of the chain contour length. 'jTgff = 0.01 
and N = 10, roughly corresponding to Wi ~ 5 — 10. All three cases correspond to the 
same far-field behavior. Finite concentrations at Z = in (a) and (b) occur due to a lack 
of divergence in equation [4] and the neglect of a short-range excluded- volume potential at 
the surface, (c) Diagram demonstrating roughly where the far-field approximation is valid 
and where the near-surface effects need to be considered (extended-chain theory) . 
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for near-surface hydrodynamics. 

We can also consider flows in a slit, where two walls are present. If the slit height H 
is on the same order of magnitude as the polymer contour length Na, we demonstrate 
extremely large differences between existing theories and our elongated chain model. Since 
we have introduced an expression for Fl{Z, N, Ft) that is based off a single image chain, 
we simply need to consider an infinite series of image chains. For a slit of height H, this 
becomes FL,sHt = ET=o^L{jH/a + Z,N,Ft) - FlUJ + l)H/a - Z,N,Ft). We plot a 
sample distribution function in a slit for our chain elongation model and both the far- 
field and non-far-field dumbbell models in Figure [Sja, where the far-field behavior is held 
constant. While the two dumbbell cases are essentially identical in this regime, there is a 
marked difference in the distribution function that results from considering the elongated 
chain geometry. 

These depletion layer calculations emphasize where our elongated-chain theory provides 
a meaningful improvement over dumbbell-based models, which is at high-shear and near- 
surface conditions. Higher shear rates result in a longer effective chain, which increases 
the height Z* where there is a crossover from the near-surface to the far- field regimes. 
This is demonstrated schematically in Figure [3|;, which indicates the extent of the near- 
surface regime as a function of flow rate (shown here in terms of the Weissenberg number 
Wi = 'jTz, for ease of comparison with traditional literature). 

In summary, we have developed an analytical expression for the lift force on a poly- 
mer near a surface that significantly improves upon existing theories by accounting for 
extended geometries. Computer simulations incorporating the hydrodynamic forces were 
performed and confirmed these analytical expressions, which give rise to the appearance 
of non-monotonic lift force behavior. This non-monotonic lift force behavior appears in a 
near-surface, highly extended regime, and it is not captured by a dumbbell theory. Further- 
more, our results have important implications for slit flows in microfluidics and polymer 
adsorption, where the distance of the polymer from the surface and its contour length are 
on the same order of magnitude. piH] Figure [3|) suggests that this result could be verified ex- 
perimentally using the visualization of the cross-channel distribution fluorescently-labelled 
polymers by applying strong shear flows in narrow channels where this effect is particu- 
larly large. Cross slit flows may also provide a tool to directly examine the lift-force on 
individual polymers. 

We acknowledge funding from the National Defense Science and Engineering Graduate 
Fellowship 

1 Supplemental Information 

1.1 Hydrodynamic Mobility Tensors 

To describe the effect of a force applied to a bead at rj on the surrounding fluid flow at r j , 
we define a mobility tensor Hij{ri,rj) that satisfies the equation Vj(rj) = fiij{ri,rj) ■ Fj(rj). 
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There are multiple forms for this tensor, depending on the assumptions made and the 
system of interest. The simplest version, called the Oseen tensor, describes the effect of a 
single point force on the surrounding fluid flow in an infinite medium: [1] 

/^ij,o(ri-rj =r) = — ^[I + ff] (5) 
Sirrjr 

where we define r as the magnitude of r and f = r/r. To consider beads of finite size, we 
can instead use the multipole expansion of the Oseen tensor, known as the Rotne-Prager- 
Yamakawa tensor: p!6l II 3| 
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where /io = l/(67rr/a) is the Stoke's mobility of a bead of radius a. 

Both the Oseen and Rotne-Prager-Yamakawa tensors assume an infinite medium. If 
there is a surface present, additional terms need to be added to these tensors such that 
they account for the no-slip boundary condition. The resulting tensor is called the Blake 
tensor, and is given (for the Oseen-based Blake tensor) by:[lj 



where rj^i^z 



+2rlzfJ'ij,PD{r + Tj^z) - 2rj^zfJ'ij,SDir + rj^z) 
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and 



^^ij,SD{r) = rzliij,PD{r) + (1 - '^^x„z) 
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1.2 Simulation Methods 

To simulate a single chain in the geometry specified in the paper, we use Brownian Dy- 
namics (BD) simulations. Our polymer is represented by a bead-spring model, which is 
composed of 2N beads i at positions rj that are held together by a harmonic potential: 
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The beads interact with other monomers through a Lennard-Jones potential: 



2N 



12 



(10) 



(11) 



10 



There is also a harmonic potential that fixes the polymer at a given height Z from the 
surface: 

Uz = '^'l^ ' (12) 

where all values designated with a tilde are dimensionless, with distances normalized by 
the bead radius a, energies normalized by fcT, and times normalized by the characteristic 
diffusion time r = Gnrja'^ /{kT). fij is the distance between beads i and j, 2N is the overall 
number of beads in the chain, k = 500 is the bead-bead spring constant, and e is a bead 
interaction parameter that controls the strength of bead-bead attraction. For this paper 
we use the value e = 0.41, which is typical for a 9- polymer. Beads move through this 
potential via integration of the Langevin equation: 

-.fi = voo(fi) - {hj,Byvfitot{i) + Vrj • D»j) + (13) 

3 

where Voo(fi) is the undisturbed solvent flow profile, /Ujj is the Rotne-Prager-Blake mobility 
matrix, Utot = Us + Ulj + Uz, Djj = ksTfiij is the diffusion tensor, and S,i is a random 
velocity that satisfies (^j) = and {Cj{t)^i{t')) = 2kBTfiij6{t — t'). Two cases are considered 
in this letter for Voo(fi). We represent the force loading corresponding to a shear flow with 
Voo,x(fi) = 7''"(ri,x — X^j'^ (2A'')), which captures the linear flow increase away from the 
center of mass of the chain. The pulling case is given by Voo^x(fi) = FQ5{i — 2N) — FQ6{i), 
which essentially only places a force in the x-direction on the flrst and last beads in a chain. 
The Langevin equation is discretized by a time step of At = 10"'^, and for a given set of 
conditions 10'' simulation steps are used. For a given value of Z, the average lift force Fl 
on a polymer is given by: 

Erfjz/(2iv)-^ 

Fl = , (14) 

5 

1.3 Polymer Dimensions in Shear Flow 

The landmark work by DeGennes on chain behavior in strong flow flelds allows us to 
calculate chain dimensions as a function of shear rate, We desire a relation that yields an 
approximate " aspect ratio" of a polymer chain in shear flow that serves as the basis for the 
respective number of beads that would be used in an analogous bead-spring representation 
of the chain. 

We begin by introducing a series of equations that DeGennes derives from the Peterlin 
Dumbell formalism: [3] 



1 - EC^^ + jTzC^ = (15) 
-EC^, + ^^TzC = (16) 

ECzZ = ECyy =1 (17) 
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where E is a chain extension parameter that is 1 at low extensions and diverges at high 
extensions, Cij = 3{rirj) /Na^ is a symmetric matrix describing the polymer dimensions, 
and 'jTz is the dimensionless shear rate (also known as the Weissenberg number) based on 
the Zimm relaxation time tz- Equations [15] and [16] can be reorganized to give the result: 
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The aspect ratio of the chain geometry, which will become 2N ^ can be then obtained using 
the relationship: 
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(19) 



This equation yields intuitive behavior - at jrz < 1, the polymer is roughly spherical. As 
jTz becomes greater than 1, there becomes a regime of linear scaling of the aspect ratio 
with the shear rate. This regime flattens out when E > 1 at high extension, and the value 
of 2A'' asymptotically approaches the contour length of the polymer. 



1.4 Derivation of Lift Force for an Elongated Chain (Pulling Scenario) 

In this letter we have derived the behavior using the kernel (dFT/dn) ~ Gnija^'y sin {26)n, 
however others could be applicable. For example, if there is an abundance of mass at the 
chain ends, it may be more appropriate to use a pulling profile (dFT/dn) = FQ5{n — N). 
This is essentially the dumbbell result, only we now integrate over the entire chain contour 
in the same fashion as before: 



Performing the integration, we get the result: 
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This simplifies to the result as Z goes to zero: 



and as Z goes to infinity: 
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Figure 4: Graph of Fl versus Z for a shear load profile. Solid lines show theoretical re- 
sults, dashed lines show the effect of bead-discretization at low values of Z, and dotted lines 
demonstrate the far-field dumbbell results. Simulation data is also shown, with filled sym- 
bols representing data without fluctuations and open symbols representing data including 
fluctuations. These results collapse onto a single curve (inset) using the scaling Z — )• Z/N. 



We can still define a characteristic height Z* where there is a transition between these two 
behaviors. Solving for the maximum of the shear behavior, we obtain the result Z* ~ A^. 
Rearrangement of the result for the pulling case can be plotted as F versus Z/N. Results 
for this are shown in Figure 1 (SM), which is analogous to Figure 2 in the letter. 



1.5 Low-Z Regime 

The replacement of the summation in equation 2 with an integral is valid if the chain 
is continuous or if the chain is in the far-field limit. If the chain is discrete and close 
to the wall, like in our simulations, this approximation does not hold. If we retain the 
summation in the low-Z limit, we can match our simulation results and demonstrate that 
this is the reason for the disparity between simulation and theory in Figures 2a and 2b. 
The discretized version of equation 2 yields the result: 



Fl. 



9Z3 JX. 



disc 



16N ^ 

n.m=0 



(n — m) 



Ftu 



[n — m)^ + 



5/2 



+ 



(n + m) 



(n + m)^ + Z"^ 



5/2 



(24) 



We can simplify this for small Z by considering that the effect of a given dipole on its 
surroundings is both antisymmetric and quickly approaches zero. Therefore, beads in the 
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center of the chain will not exert a strong net lift force. Only the beads close to the chain 
ends contribute non-negligibly to the lift force, which we express as a truncated series: 



^ L ,disc,shear ^ ^ / ^ 

i=0 
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(25) 



for shear and: 
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for pulling. I is a small number representing the point at which the summation is truncated 
(we use / = 3). This result agrees with the values obtained for simulation data at small Z, 
as shown by the dotted lines in Figures 2 and 1(SM). This indicates that the main reason 
for the disparity between theory and simulation is due to the discretization of the chain 
in the simulations, so therefore we expect that a more realistic chain would more closely 
agree with the original, continuous theory. 



1.6 Translating From Geometric to Chain Parameters 

The theory presented as such is entirely based on geometric parameters, such as the polymer 
aspect ratio N and the vertical length scale a. In this respect, no particular chain model is 
assumed and the theory is very general. It is possible to incorporate chain characteristics 
into such a theory through a change of variables. We take the De Gennes model of a 
polymer in flow as a simple example, and we use scaling to show how we can obtain 
the traditional scaling relationship found for the far-field dumbbell model. [3] We use the 



far- field result for shear flow, F^^shear ~ §7''" sin (26*)^, and make the following scaling 
replacements: 



N r^Wi 
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sin (20) ~ Wr^ 



(27) 
(28) 
(29) 



which are appropriate at intermediate values of Wi. [3] n. is the degree of polymerization, 
and b is the monomer size. The translation of r to tz is done by noting that r ~ "qo? / (kT) 
and Tz ~ rfJT'rv'/'^ / {kT). This leads to the result: 



Fl 



(30) 
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where all distances are now scaled by b. The depletion zone can be calculated by coun- 
terbalancing the lift force with the diffusive force pushing the polymer back towards the 
surface, Fs- Since there is a no- flux condition at the surface, we expect the concentration 
profile to be of the form c{Z) = c(0) + AZ^ where A is an arbitrary constant. The diffusive 
force is thus given by Fs ~ Vc(Z) ~ —Z. There is then a critical height Z* where these 
forces balance, which is the depletion length: 

Z* ~ 1^^2/3^1/2 (31) 

which is the relationship given in the literature. [15j We add the caveat that these scaling 
relationships are only appropriate in a finite range of Weissenberg numbers, and is thus not 
universal. Separate analysis would have to be performed to obtain this sort of relationship 
under different conditions, however the underlying geometric theory is completely general. 

In the letter, order of magnitude calculations are done to relate the geometric param- 
eters to the chain parameters. This is done using the relationship (from the above scaling 
analysis) : 

— z ^ 

Using the values and relationships = 10, 'JTefj = 0.01, and n^^'^b/a ~ Wi^^^, we obtain 
the result: 

Wir^b (33) 

This is an order of magnitude result, but it demonstrates that the results presented in this 
paper are in the relevant shear rate regime to see the effects indicated (lift and depletion). 
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